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An explanation of the dynamical mechanism for apse- 
alignment of the eccentric Uranian rings is necessary before 
observations can be used to determine properties such as 
ring masses, particle sizes, and elasticities. The leading 
model (Goldreich and Tremaine 1979) relies on the ring 
self-gravity to accomplish this task, yet it yields equilib- 
rium masses which are not in accord with Voyager radio 
measurements. We explore possible solutions such that the 
self-gravity and the collisional terms are both involved in 
the process of apse- alignment. We consider limits that cor- 
respond to a hot and a cold ring, and show that pressure 
terms may play a significant role in the equilibrium condi- 
tions for the narrow Uranian rings. In the cold ring case, 
where the scale height of the ring near periapse is compa- 
rable to the ring particle size, we introduce a new pressure 
correction pertaining to a region of the ring where the 
particles are locked in their relative positions and jammed 
against their neighbors, and the velocity dispersion is so 
low that the collisions are nearly elastic. In this case, we 
find a solution such that the ring self-gravity maintains 
apse-alignment against both differential precession (m = 1 
mode) and the fluid pressure. We apply this model to the 
Uranian a ring, and show that, compared to the previous 
self-gravity model, the mass estimate for this ring increases 
by an order of magnitude. In the case of a hot ring, where 
the scale height can reach a value as much as fifty times 
larger than a particle size, we find velocity dispersion pro- 
files that result in pressure forces which act in such a way 
as to alter the ring equilibrium conditions, again leading 
to a ring mass increase of an order of magnitude; however, 
such a velocity dispersion profile would require a different 
mechanism than is currently envisioned for establishing 
heating/cooling balance in a finite-sized, inelastic particle 
ring. Finally, we introduce an important correction to the 
model of Chiang and Goldreich (2000). These authors re- 
lied on collisional forces in the last ~ 100 m of a ~ 10 km 
wide ring to increase ring equilibrium masses by up to a 
factor of ~ 100. We show that their analysis leads to a 
strong, artificial ring mass dependence on the adjustable 
parameter A (the lengthscale over which the ring’s optical 
depth drops from order unity to zero at the edge). A cor- 
rected treatment of ring edges that takes into account their 
ridge-like structure retains the increase of ring mass of the 
order of ~ 100 for a 10 km wide ring, while exhibiting 
weaker A dependence. We introduce shepherd satellites 
and show that they can have an effect on equilibrium ring 
masses and surface density profiles. We conclude it is likely 
that a modified CG model can account for the masses of 
narrow, eccentric planetary rings; however, the role of 
shepherd satellites both in forming ring edges and in alter- 
ing the streamline precession conditions near them needs 
to be explored further. It is unclear whether such a model 
allows for the possibility of rings with negative eccentricity 


gradients. 


1 Introduction 

The nine main narrow Uranian rings are each less than 100 
km wide, with moderate optical depths (> 0.3), and mostly 
made of particles larger than centimeters. Most have been 
found to be eccentric, inclined, and sharp edged (French et al. 
1991). Because of Uranus’ oblateness (J 2 = 3.343 x 10 -3 ) 
the pericenter of a ring particle will precess at a rate dependent 
on its semimajor axis. As a result, particle orbits at the inner 
and outer edges of the narrow ring should quickly misalign. 
Therefore, the presence of such eccentric rings speaks unam- 
biguously to the need for a mechanism to counter differential 
precession. 

The leading theory to explain ring equilibrium uses the 
ring’s own self-gravity to enforce ring particle apsidal align- 
ment (Goldreich and Tremaine 1979). However, this theory 
has run into trouble when confronted with observations, and 
has failed to predict the dynamical behavior of ring streamlines 
or even to explain the overall state of some rings (Graps et al. 
1995; Goldreich and Porco 1987). In particular, the surface 
density implied by the self-gravity model is smaller than seems 
plausible for the a and /3 Uranian rings based on the apparent 
particle size (Tyler et al 1986; Gresh 1990). Furthermore, it 
has been pointed out that with the low surface densities im- 
plied by the self gravity model the torques exerted by inner 
shepherd satellites would be insufficiently strong to confine 
these two rings against drag from the distended exosphere of 
Uranus (Goldreich and Porco 1 987). Unfortunately, alternative 
apse-alignment models have not been developed quantitatively 
(Dermott and Murray 1980), or have been based on unrealistic 
assumptions about the viscous stress tensor (Papaloizou and 
Lin 1988). 

Recent work by Chiang and Goldreich (2000) proposed 
that the collisional forces near ring edges can increase the 
mass of a 10 km wide ring by up to a factor of ~ 100, making 
it possible to reconcile the observational mass estimate for this 
ring based on radio occultation constraints of ring particle sizes 
with the dynamical mass estimate based on the lack of differ- 
ential precession. Their model allows for both positive and 
negative ring eccentricity gradient, which is not easy to recon- 
cile with the the observation that all narrow rings have positive 
eccentricity gradients (considered as evidence in support of the 
self-gravity model in the first instance). 

In section 2 we look at the case of a cold ring in which 
stresses build up due to particle jams (regions of the rings where 
particles are locked in their relative positions with neighboring 
particles). In section 3 we increase the velocity dispersion of 
the ring and set constraints on its azimuthal behavior assuming 
that collisional stresses alter the equilbrium conditions for apse 
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alignment. In section 4 we modify the model of Chiang and 
Goldreich (2000). In section 5 we present our conclusions and 
discussion. 


2 Cold Ring 

A cold ring is one in which the scale height of the ring near 
periapse is comparable to the ring particle size. Here we model 
the effects of a ring particle jam on the overall dynamics of the 
ring. We define a “jam” as a region of the ring where particles 
stay in contact with each other, such that it is not valid to treat 
particle collisions as pair-wise hard sphere interactions. As we 
shall see, such a region would be characterized by enhanced 
pressure terms. 

Other publications have advanced the possibility of such 
a state (Wisdom and Tremaine 1988). Unfortunately, they 
contain unrealistic assumptions about the ring, such as the 
treatment of ring particles as smooth spheres. The conditions 
under which more realistic ring particles jam are not well 
understood. Nevertheless, the increase in ring self-gravity 
and optical depth as the ring converges will lead to a rapid 
increase in the number of inelastic collisions and lower the 
velocity dispersion (Mosquetra 1996), thus creating conditions 
which may lead to a ring particle jam. For the puipose of this 
paper, we will assume that a particle jam can take place in a 
converging flow, and we will simply parameterize its vertical 
extent in terms of the local ring optical depth. 

2.1 Fluid Description 

Streamlines in a perturbed ring can be characterized by an 
azimuthal wavenumber of the distortion m > 0 and a pattern 
speed ftp. In a frame of reference which rotates with the 
pattern speed we have r = a(l — e(a)cos0), where r is 
the cylindrical radius, a is the semi-major axis, e << 1 is 
the eccentricity, 0 = m(<j> + A), and A is the phase angle. 
The Jacobian of the transformation between (r, 0) and (a, 9) 
is given by J(a,8) = 1 - geos (0 + 7 ), in general q e — 
geos 7 = a^, and g 7 = q sin 7 — moe^. From here on 
we will assume that 7 < < 1 and q ~ g e . 

In our model, the pressure in the jammed ring region is 
not the result of distinct collisions. In the limit as collisions 
become very numerous and gentle, we resort to a pressure 
description wherein a given ring particle spends most of its 
time in contact with its neighbors. We find the following 
expressions for the pressure and viscosity in this tightly packed 
region with low velocity dispersion a 

po~ 0.1 ( 1 ) 

7,0^ 0.01 k c E*(^f-y, ( 2 ) 


where R and p p are the particle radius and density, k c is the 
viscosity flux parameter which measures the relative strength 
of the fluid pressure and viscosity, and E ~ 10 11 dyne cm™ 2 
is the elastic modulus of ice at temperatures comparable to 
those of ring particles (Borderies et al. 1984) (in this context 
particle collisions produce nearly elastic deformations of ice 
(Hertz 1881)). Physically, these equations result from the 
redistribution of the stress taking place at the area of particle 
contact over their surface area. 

Furthermore, the vertical pressure does not simply balance 
the vertical acceleration of the fluid (Borderies et al. 1985). 
If that were the case the ring would simply puff up, and the 
stresses would never reach the values given by Eqs. ( 1 ) and (2). 
Instead, we assume that, so long as ring stresses do not exceed 
the yield modulus of ice, ring particles can be prevented from 
moving to a higher plane by the surface roughness of their 
nearest neighbors. That is, the jam itself can prevent particles 
from moving not only in the ring plane but in and out of it as 
well. By contrast, in a smooth sphere model the difference in 
the vertical location of particle neighbors inevitably results in 
forces that push particles to higher vertical positions, which 
acts to hinder the formation of jams. 

Next, we find the height integrated pressure tensor com- 
ponents 

P rr = 2 h (p 0 ~ , (3) 


„ ^1 — 4/3gcos0 

PtB — 3 / 1770 ^ j (4) 

where Q = (GM p /a s ) 1/2 is the orbital frequency, M p is the 
planet’s mass, a is the semi-major axis, and h is the jam’s 
scale height. The optical depth of such a ring is given by 
r = ro/J, where ro is the optical depth at quadrature. As 
the ring converges, particles get closer together and the ring 
density increases. Around unit optical depth it is possible for 
particles near the midplane to jam. If the contraction continues 
from that point a larger fraction of the ring, of vertical extent h, 
will become jammed. We model this process by the heuristic 
equation 

where R is the particle radius. Roughly speaking, the above 
equation simply describes the vertical extent or thickness of 
the jammed region. Since not all ring particles need be in 
the gridlocked vertical portion of the ring, this scale height 
should not be confused with the actual ring scale heighj H y 
which is determined by the combination of the ring’s velocity 
dispersion and the exclusion of volumes in the ring (Mosqueira 
1996). The density of the ring within the jammed region 
h is assumed to be the maximum attainable ring density or 
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p = p P FF ~ 0.33 g cm -3 , where p v ~ 0.5 g cm -3 is 
the density of a ring particle and FF is the maximum ring 
filling fraction. It should be noted that, just as in the case 
of the traffic jams we are all used to, the ring particle jam 
will be over near 0 = 0, when the divergence of the mean 
flow turns positive. That is, for simplicity we assume that the 
jam takes place between the negative angle 0o and 0, where 
cos 0o = (1 — To)/q and (1 — q) < to < 1. 

Next, we use a two- streamline model of the ring (Borderies 
et ai 1983) and define the viscous stress frequencies 
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where Aa is half the thickness of the ring, and M r is the mass 
of the ring. Finally, the equilibrium conditions for the ring are 
given by (Borderies et ai 1983) 
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where 2Azu° is the ring’s apsidal shift, q = aAe°/Aa, and 
the libration frequency due to the ring’s self-gravity is 
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Figure 1: Apsidal shift in degrees as a function of mass 
for the a ring that satisfy the equilibrium conditions (eqs. 
(7) and (8)) for several choices of the viscosity coefficient 
k Cy which determines the relative strength of the pressure 
and the viscosity. For small k Cy larger masses correspond to 
a more extended “jam” in the ring. As the jam lengthens 
beyond 0 = 0 (the solution corresponding to larger values 
of q) } the P rr term begins to dominate the calculation of Ai 
and this leads to additional solutions at higher masses. For 
the choice of k c = 1.2, a solution is found with a negative 
apsidal shift (~ — 1°) and a mass (~ 20 M sg ) that yields 
a surface density consistent with the value obtained using 
S occ = | ppRr ~ 30 g cm -2 with R = 100 cm. For the a 
ring e = 7.81 x 10~ 4 , q = 0.45, and tq = 0.55. 


where 


H(q 2 ) 


1 - y/1 - 
q 2 yf\ - q 2 


( 10 ) 


For the case that n., ff >> |Ai|, |A 2 | (negligible pressure ef- 
fects) we obtain the self-gravity mass for the a ring M sg = 
3.4 x 10 16 g 1 using a perturbation parameter for the ring of 
q ~ 0.45 and eccentricity of e = 7.81 x 10~ 4 . This corre- 
sponds to a surface density ~ 2 g cm' 2 , which is about an 
order of magnitude smaller than the surface density given by 
Soec = 4/3 p p Rr ~ 30 g cm 2 , using R = 100 cm and 
ro = 0.55, as favored by radio occultation observations of the 
Uranian rings (Tyler et al. 1986; Gresh 1990), and p p — 0.5 
g cm -3 . Note that we implicitly assume that the elastic prop- 
erties of the a ring particles are those of solid water ice. 

^ote that this value for the ring mass was obtained using a 2 
streamline model. As we will see later in section 4, the use of more 
streamlines yields a smaller mass 


In the case of the a ring, the lack of any measurable 
differential opacity between wavelengths of 3.6 cm and 13 
cm indicates a relative depletion of centimeter sized particles 
(Gresh 1990). Furthermore, the proximity of this ring to the 
planet and the extended nature of the planetary atmosphere 
would make sub-centimeter sized particles short-lived (Gol- 
dreich and Tremaine 1979). These arguments make it appear 
very unlikely that the surface density of the a ring could turn 
out to be as low as the self-gravity model indicates. This cir- 
cumstance leads us to search for larger mass solutions to the 
equilibrium Eqs. (7) and (8). In particular, we are interested 
in solutions such that A 2 >> |Ai| > 0 and Ai < 0. Such 
solutions would make it possible to increase the dynamical 
mass estimate of the ring while keeping the streamlines nearly 
aligned (with a small negative apsidal shift as observed). Our 
new pressure terms allow for them. 

To find the contribution of the “jammed” region to Ai and 
A 2 we estimate the velocity dispersion a in the region between 
the negative angle 0o and 0, where the jam occurs. According 
to Andrews’ theory (Andrews 1930), collisions are nearly elas- 
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tic when the impact velocity is v < v* ~ 0.03 cm s“ 1 . This is 
also the approximate impact velocity at which the yield stress 
of ice is reached (Borderies et ai 1984). Numerical experi- 
ments indeed point to very low impact velocities in a ring (Wis- 
dom and Tremaine 1988; Richardson 1994; Mosqueira 1996; 
Salo 1991) given a velocity dependent coefficient of restitution 
e (Bridges et al. 1984), which is the factor by which the rela- 
tive normal velocity is reduced after a collision. This velocity 
turns out to be of the order of a large particle’s escape velocity. 
For a particle with radius R = 100 cm and density p p = 0.5 
gm cm -3 this velocity is ~ 0.05 cm s” 1 . For the a ring a 
typical relative velocity between two neighboring particles of 
this size due to the velocity gradients in the fluid flow is of 
order — 0.02 cm s” 1 . This means that both the ring shear and 
its convergence need to be considered when calculating the ve- 
locity dispersion a (Brahic 1977). Because the jam takes place 
in the ring region where the material is converging, we used 
a 1 = v+ + (RQq sin 0/ J) 2 4-/? 2 n 2 (16g 2 — 15+24J)/(8J 2 ). 


2.2 Cold Ring Results 

Since the jam takes place just before periapse where 9 < 0, the 
calculation for A 2 involves a positive integrand P rT cos 0, thus 
Eq. (6) shows that the jam’s contribution to A 2 is positive. In 
general, A 2 need not be small when compared to the libration 
frequency Vt ag . This means that the jam can significantly alter 
the calculation of the mass in the self-gravity model, contrary 
to earlier treatments of the pressure terms. In fact, a positive 
A 2 means that a given ring perturbation will require a larger 
ring mass to satisfy the equilibrium condition given by Eq. 
(8) The case of Ai is somewhat more complicated. For small 
negative angles the calculation of Ai is determined by the 
positive P T o cos 9 term, but as the jam increases in extent it 
becomes possible for the P rr sin 9 term to provide a negative 
contribution that dominates the calculation of Ai (hence the 
turnover in the small k c curves in Fig. 1). It is this behavior of 
Ai which makes the desired solution A 2 >> | Ai | with Ai < 0 
possible. 

Because the jam contribution to A 1 is not necessarily large, 
we need to consider the effects of the rest (non-jammed por- 
tion) of the ring. That is, we need an estimate of the apsidal 
shift provided by the ring as a whole. Such a calculation is 
model dependent, but the choice does not affect the thrust of 
the argument presented here. We choose an apsidal shift in the 
absence of a particle jam of ~ —3.5°. We simply add the Ai 
corresponding to this negative apsidal shift to the "jam” contri- 
bution, and we find the solutions of the equilibrium equations. 
These solutions turn out to be rather sensitive to the (unknown) 
relative strength between the fluid pressure and viscosity ( k c 
in Eq. (2)). We choose some representative parameters and 
show that it is possible for the actual ring mass to be many 
times that of the self-gravity model. It should be noted that 
in all fitted cases observational results yield negative apsidal 
shifts which can be used to constrain the value of k c . 

For small k c , larger masses correspond to a more extended 
“jam” in the ring. As the jam lengthens for 9 < 0 away from 


periapse (the solution corresponding to larger values of q ), the 
P rr term begins to dominate the calculation of Ai and this 
leads to additional solutions at higher masses. For the choice 
of = 1.2, a solution is found with a negative apsidal shift 
(~ -1°) and a mass (~ 20 M sg ) that yields a surface density 
roughly consistent with £ 0 cc, which is more than an order of 
magnitude larger than = M sg /(4naAa) ~ 2 g cm^ 2 . 
The perturbation parameter q for the ring fell between 0.45 
and 0.5, which corresponds to a maximum h of 0.1 R. The 
jammed ring region occurred between 0° and a minimum of 
- -25°. 
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Figure 2: Apsidal shift as a function of mass for the a ring 
for several larger values of the viscosity coefficient. Here we 
demonstrate that it is possible to relax our assumptions on k c 
and obtain lower mass solutions (still several times the self- 
gravity mass) for the ring. In order to obtain these solutions, 
we must choose an alternate form of the scale height (Eq. 
(11)) where the peak stress values occur away from closest 
approach. This essentially means that the largest fraction of 
particles that become jammed occurs far from 0 = 0. Larger 
masses correspond to a more extended pinch in the ring, and 
larger values of q . Masses larger than roughly 6 times the self 
gravity mass (long-dashed curve) are not possible with this 
form of the scale height since the integrals over the pressure 
tensor approach a limiting value. As one increases k c , h 
becomes more sharply peaked. 

2.3 Other Jam Models 

In order to explore the model dependence on the specific form 
of the heuristic equation Eq.(5), we use the expression 

h = R[J-l]{e/Bo) n , ( 11 ) 
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where n is a positive integer. The idea behind this formula 
is to shift the peak position of largest stress. By doing this, 
the model allows for larger turnover values of k c while still 
significantly increasing the ring’s mass. In Fig (2), we plot the 
corresponding results using this new heuristic expression. It 
should be noted, however, that the mass increase is consider- 
ably less than in our earlier model. 


3 Hot Ring 

Here we consider the possibility that the increased collisional 
stresses result from larger velocity dispersions in the fluid than 
result from a straightforward balance of viscous heating terms 
and inelastic collision cooling terms (Goldreich and Tremaine 
1978). So long as the Bridges et ai coefficient of restitution is 
assumed, equilibrium velocities on the order of ~ 0.1 cm s _1 
result, which are too small to lead to stresses that can affect ring 
precession. Nonetheless, there are mechanisms that can con- 
ceivably change the energy balance for a narrow ring such as 
the presence of shepherd satellites, the presence of ring struc- 
ture, large particles, and even global effects (Shu et al. 1985). 
Here we search for a profile of the velocity dispersion that can 
produce stresses leading to an order of magnitude increase in 
the ring mass estimated for the a ring. We use the pressure 
and viscosity expressions for local transport (Mosqueira 1996) 


po = 0.33pcr 2 , (12) 


770 = 0.033p p cri?, (13) 


where p = E/(2 H) is the ring density, and we obtain the scale 
height for the ring from the equation for hydrostatic balance 


H = (-ttGE + ((ttGE) 2 + ft V/ 3) 1/2 ) /ft 2 , (14) 


which takes into account the ring self-gravity (Mosqueira 
1996). 
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Figure 3: Velocity dispersion profiles for the a ring. These 
represent two examples of profiles that yield values for the 
viscous frequencies Ai and A 2 such that the surface den- 
sity of E = 36 g cm' 2 , which is consistent with E occ . In 
the solid curve ( a oc 1/J), the peak velocity occurs just 
before 0 = 0. The dispersion velocity varies over the en- 
tire ring with the minimum velocity occuring near apoapse. 
In contrast, the dotted curve demonstrates that a velocity 
dispersion with the opposite behavior may also fit the ring 
mass. Here the maximum velocity occurs after 0 = 0 and 
the minimum occurs just before periapse. However this fit re- 
quires relatively little variation in velocity far from the apse, 
and a steep transition from close densely packed particles 
with dampened random velocities to the maximum disper- 
sion velocities after 6 = 0. The resulting apsidal shift is 
AtJ7° ~ —2° and the ring mass is ~ 20 M sg . 

3.1 Hot Ring Results 

We find that, for sufficiently large velocity dispersions a, the 
stress tensor can change the equilibrium conditions Eqs. (7- 
8) in such a way as to lead to a ring mass consistent with 
a surface density of ~ 40 g cm -2 for the uranian a ring, 
while still preserving ring alignment. In Fig. (3), we plot 
two velocity dispersion profiles for the a ring that satisfy the 
above requirements. Namely, both of these profiles yield the 
same values for Ai = —7.63 x 10“ 10 s~ l and A 2 = 6.30 x 
10" 8 s _1 . Furthermore, A 2 ~ = 6.52 x 10' 8 s“ l . 

While the dashed line satisfies the required equilibrium 
conditions, it is a very unlikely profile and we will have nothing 
further to say about it. For the solid line, the peak occurs before 
periapse, and the scale height (Eq. (14)) can be as high as 50 
times the particle size. This velocity profile is difficult to 
justify on the basis of simple heating/cooling models. The 
high velocity dispersion in regions of high optical depth are 
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difficult to reconcile with cooling due collisions, since regions 
of high optical depth lead to many more collisions. We tried 
the presence of large particles as a heating mechanism that 
would lead to extra heating in high optical depth regions, but 
found that large particles are unlikely to produce the required 
strong contrast between the velocity dispersion of high and 
low density regions. There are, however, mechanisms, such as 
heating due to ring structures forming in high density regions, 
which might lead to such a velocity dispersion, but they are 
difficult to characterize and we will not attempt to do so here. 


leads to ridge-like edges with strong density gradients present 
on both sides of the density peak at each ring edge. As a 
result, the surface density profiles in CG are inconsistent with 
the assumption that went into producing them. Unlike CG, we 
do not assume a linear collision term for the edge streamlines, 
instead we use a boundary condition that fixes the relative 
masses of edge streamlines and include a collision term based 
on the local surface density gradient for all other streamlines. 
The collisional acceleration for a ring particle is given by 


4 Edge Effects 


VP al rfS . 

£ (1 — q e cos /)£ da r> 


(17) 


Recently, Chiang and Goldreich (2000) (hereafter CG) pointed 
out that collisional forces felt by material in the last ^ 100 m 
of a ~ 100 km ring can increase the equilibrium masses up to 
a factor of ~ 100. If correct, their study can serve to reconcile 
the observational constraints on ring surface densities with 
theory, though it would raise anew the issue of why all narrow 
planetary rings exhibit a positive eccentricity gradient. 

Here we take a look at their model and show that their 
model treats ridge-like edges improperly. Like these authors, 
we divide the ring into N equally spaced, apse-aligned stream- 
lines with constant q e = aAe/Aa. The j th streamline has 
mass mjy semi-major axis aj = « + — (N 4- l)/2] Aa/N , 

and eccentricity tj = e + \j — ( N 4- 1 )/2] A e/N. For the 
ring to maintain uniform precession the effects of the planetary 
oblateness, self-gravity and interparticle collisions must coun- 
teract one another. The contribution to differential precession 
due to the planetary oblateness is 



that of self-gravity is given by 


A; 



qe-tf(g g 2 ) f , a E m k 
7re A/p 3 aj — ’ 


(16) 


where H(q 2 ) is found from eq. (10). 

Also like CG, we ignore the decrease in velocity dispersion 
from the ring edge cr*, ~ 1 cm s" 1 to the interior. This 
assumption helped our numerical solutions to converge and 
it is justified because the width of a satellite resonance is 
of order w r ~ a(M a /M p ) 1 ^ 2 ~ 1 km, which covers the 
region of the ring where surface density gradients are strong 
(though some of our results would be marginally altered by 
using a lower velocity dispersion Ci ~ 0.1 cm s -1 in the 
ring interior; in particular, the ring surface density away from 
the edges would increase). CG treated each ring edge by 
introducing a linear collision term for edge streamlines. Their 
treatment, however, implicitly assumes that the edges are one 
sided surface density drops, whereas the result of their analysis 


where P is the height-integrated pressure and / is the true 
anomaly. We then have 


(dw\ qJHfool 1 dZ 

\ dt / c Qae £ da 

The collisional term is then given by 


(18) 
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if \ < x = (j — 1 / 2 ) Aa < 4*. 

if\<y = (N-j + l/2)fr < Aa; 

(19) 


Given the mass of the streamline just outside each bound- 
ary we fix the mass of boundary streamlines (streamlines 
within a distance A from either edge) using a simple fit given 
by mj — mj±ie~ Aa / NX t where the plus sign corresponds 
to the inner edge and the minus sign corresponds to the outer 
edge. Boundary streamlines do not enter in the system of 
equations other than as a mass term in the self-gravity piece 
and as a collision term for the two streamlines just outside each 
boundary. Thus, we treat each edge as a density ridge. On one 
side of the ridge the density drop is picked up by the boundary 
condition, and on the other side of the ridge the density drop 
is taken into account through the collision term. 

The assumption is that the shepherd satellites alter the pre- 
cession conditions of boundary streamlines in such a way as to 
produce the needed drop off while preserving apse alignment. 
Ultimately, proper treatment of the boundaries needs to in- 
clude the effects of shepherd satellites. A fully self-consistent 
calculation, however, is left for another paper. 


4.1 Surface Density Profiles 

The system of equations to be solved in the CG treatment are 
linear in nature and can be readily solved by the use of matrix 
inversion techniques. The CG system involves N equations 
in N unknowns for the streamline masses mj which we solve 
by way of an LU decomposition of the matrix of coefficients 
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determined using the self-gravity term (Eq. 16). The decom- 
position of this matrix into “lower” and “upper” triangular 
matrices allows us to take advantage of the relative ease of 
solving a triangular set of equations by forward and backward 
substitution. In order to reproduce the results of CG for param- 
eters appropriate for the a ring, we typically used N = 3000 


streamlines to ensure convergence. 
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Figure 4: The CG solution for the a ring with <7b = 1 cm 
s _1 is used here to demonstrate the inconsistent treatment 
of the collisional precession terms near the ring edges. The 
top panel shows the a ring surface density at quadrature as 
a function of distance from the from the ring’s centerline. 
The derivative of the surface density is taken from the CG 
solution and used to generate a new surface density profile. 
The top panel shows the original CG profile (dotted line). 
The derivative of the surface density taken from the CG 
solution is used to generate a new surface density profile, 
which is then plotted (solid line). The bottom panel contains 
the collisional precession terms used to generate each surface 
density profile. The collisional terms are qualitatively and 
quantitatively different leading to different surface density 
profiles. Here q e ~ 0.47 and the number of streamlines is 
N = 3000. 

The addition of our collisional precession term to the sys- 
tem leads to a non-linear set of equations, and thus a more 
difficult problem to solve. In order to find solutions, we 
must turn to root finding techniques. To this end, we em- 
ploy the multi -dimensional secant method developed by Broy- 
den (1965) which is globally convergent in almost all cases. 
The method requires an initial guess for the mass distribu- 
tion. The actual guess need not be close to the actual solution, 
although good guesses are rewarded with much faster conver- 
gence times. This approach reproduces the linear solutions 
of CG regardless of the choice of mass distribution with little 


effort. However, some care had to be taken for the non-linear 
cases in order to avoid convergence to local minima. Once the 
general shape of the solution was known, it could be used as 
a first guess for different A, and for any choice of streamlines 
N. A tolerance parameter for convergence of 1 x 10 “ 8 was 
used to ensure good accuracy in all non-linear cases. 
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Figure 5: A comparison of results for the a ring using 
the CG treatment as well as the treatment presented in this 
study. We plot the surface density at quadrature as a func- 
tion of distance from the ring’s midline. We chose A = 500 
m for the CG solution and A = 250 m for ours. This was to 
ensure that the the density peaks were aligned. The presence 
of our collisional precession term causes the density to dip 
away from the boundaries. Note, however, that it never dips 
below the self-gravity solution which is included for compar- 
ison. We also show that both treatments allow for negative 
q e solutions (where the perturbation parameter is the same 
magnitude as the positive case). The value of q t ~ 0.47, 
Oh — 1 cm s ~\ and the number of streamlines N = 500 for 
our cases. The corresponding CG cases have N = 3000. 

To illustrate the inconsistency in the CG treatment of the 
boundary and motivate our own treatment, we take the deriva- 
tive of the surface density profile generated by their collisional 
term. Then making use of Eq. (18) we plug this derivative 
back into the matrix as a collisional term and we generate a 
new surface density profile. A self-consistent treatment would 
yield a reasonably similar profile to the original one. Fig. (4) 
shows that the CG treatment is not self-consistent in the above 
sense. The bottom panel shows the collisional term for the 
a ring (with A = cr^/H ~ 40 m) taken from CG compared 
to the derivative of the surface density profile at quadrature 
obtained using their method. The top panel shows the original 
surface density profile obtained by the CG method with the 
density profile obtained by using the derivative of surface den- 
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sity as a collisional term. It is clear that the two profiles are not 
quantitatively or qualitatively similar. Most significantly, their 
collisional term treats the ring edge as a density drop instead 
of a density ridge, as can be seen from the comparison in the 
bottom plot. 


that the CG solution is more sensitive to the specific choice 
for A that our solution. The top plot shows the reason for this. 
Increasing A makes our profiles flatter overall but the total ring 
mass remains fairly constant. This a direct consequence of our 
treatment of the collisional term. 
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Figure 6: Variation of ring mass for the a ring as a function 
of the number of streamlines N for our model. The mass 
converges fairly quickly with the number of streamlines. The 
change in mass between N — 500 nd N - 1000 is already 
at the 1% level. For this case, A = 500 m and &b = 1 cm 
s’ 1 . 

Fig* (5) is a comparison of a ring results obtained using 
CG treatment to our own and the self-gravity model. For more 
meaningful comparison we chose different A values chosen in 
such a way that the location of the density peaks align. In this 
case we used a A = 500 m for the CG treatment and A = 250 m 
for our treatment. The solid curve, which corresponds to 
our treatment, dips much lower due to the presence of the 
collisional term away from the boundaries, but it is always 
above the value provided by the self-gravity solution. Our 
treatment locates a larger fraction of the ring mass towards the 
edges than that of CG. Had we included lower values for the 
velocity dispersion in the ring interior (given by Oi ~ 0.1 cm 
s” 1 ) the shape of the solution would have changed somewhat, 
however, due to solution convergence issues we have chosen 
not to provide those results at this time. This plot also shows 
that both models allow for negative q e solutions. 

Fig. (6) shows convergence of our solution for the a ring 
with the number of streamlines. Notice the scale of the plot. 

Fig. (7) shows the results of runs where the value for 
the parameter A was changed for the a ring from a value of 
A ~ ~ 50 m to a value of A ~ 0.bw r ~ 500 m. The 

number of streamlines within the boundary was kept constant 
at ~ 10 as A was changed. In the bottom panel we can see 
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Figure 7: Variation of a ring mass with A for our model and 
that of CG. The number of boundary streamlines was kept 
constant at 10 and solutions were found for A between 50 
to 500 m. Since the number of boundary streamlines is kept 
constant, the total number of streamlines needed in each case 
is different, with a larger number of streamlines required for 
smaller A (N ranges from several hundred to over 1000 in 
the cases shown). The CG treatment is more sensitive to the 
choice of A. In contrast, our treatment leads to masses that 
are roughly constant over A. The top pane! demonstrates 
that, for our treatment, increasing lambda causes our profiles 
to be flatter overall with the overall ring mass remaining 
fairly constant. Here &b = 1 cm s" 1 . 

In Fig. (8) we place two 5 km satellites 50 km away from 
the ring with zero eccentricity. Furthermore, we average their 
positions over one orbit and we treat them like streamlines. 
The only difference between this streamline and a regular ring 
streamline is that its mass is fixed and its perturbation param- 
eter qi a = o o e j a , where a s is a satellite’s semimajor axis, 
must be computed for each ring streamline i during evaluation 
of the precession term due to self-gravity. It is worth noting 
that in this case the ring loses its symmetry. The resulting pro- 
file has a larger density peak on the outer edge. This treatment 
neglects the role of individual satellite resonances in mantain- 
ing the ring edge. A more detailed, self-consistent model is 
left for a later paper. 
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Figure 8: a ring surface density profiles for cases with and 
without satellites for our model The dotted line corresponds 
to a case with no satellites and a A = 250 m. The solid line 
is a surface density profile in which two 5 km satellites have 
been placed at a distance of 50 km from the inner and outer 
edges of the ring. Both satellites have e = 0 and are treated 
as streamlines by averaging their positions over one orbit. 
CG’s model for the same A value is shown for comparison. 
For this case, N = 500 and <Jb = 1 cm s _ . 


5 Conclusions 

In a dense, cold ring, where a ring particle jam takes place, 
the enhanced pressure between the ring streamlines (which 
occurs just before the location where their separation is the 
narrowest) acts to destroy their alignment and must be coun- 
tered by a larger self-gravity term than would be the case in 
the presence of differential precession alone, which means that 
a larger ring mass is needed if apse- alignment is to be main- 
tained. The application of our model to the Uranian a ring 
shows that, compared to the previous self-gravity model, we 
can increase the mass estimate for this ring by an order of 
magnitude. While this model is dependent on the validity of 
a poorly known parameter (fc c ), solutions can be obtained for 
a range of plausible parameters, with perhaps more realistic 
values for this parameter. It remains to be shown, however, 
whether such '’jams" occur at all. Since occultation observa- 
tions are more consistent with a constant equivalent depth than 
a constant equivalent width, they favor a many particle thick 
model (French et al. 1991). On the other hand, the jam need 
only take place over a short portion of the ring. Furthermore, 
Voyager radio observations of optical depth and phase indi- 
cate the presence of anomalous phase behavior possibly due 
to coherent interaction or close-packing (Gresh 1990). 


In the case of a hot ring, we have modeled the velocity 
dispersion as a function of azimuth, and we have shown that 
it is possible for the pressure terms in a ring with optical 
depth close to but less than unity to produce a profile such 
that the mass estimate for the ring is again increased by an 
order of magnitude. As before, the model parameters were 
selected to explain the a ring mass discrepancy. Nevertheless, 
it is not easy to identify a mechanism that could lead to the 
required azimuthal velocity dispersion profile, which has a 
maximum close to the region of largest optical depth. It is 
possible, though unlikely, that wakes or other ring structure 
developing preferentially in high density regions would pump 
the velocity dispersion there by an amount greater than the 
increased number of inelastic collisions can dissipate, resulting 
in such a profile. 

We modify the model of CG and show that proper treat- 
ment of the pressure gradient terms near ridge-like edges leads 
to equilibrium masses which are weakly dependent on the un- 
known value of A (the length scale over which the ring’s optical 
depth drops from order unity to zero). Compared to CG, our 
treatment allocates a larger fraction of the total ring mass to- 
wards the edges. The resulting equilibrium ring masses for a 
10 km ring are increased by a factor of ~ 100 with respect to 
the self-gravity model so long as the ring velocity dispersion 
at the ring edge is about ~ 1 cm s“ 1 . Because heating/cooling 
collisional balance is likely to lead to significantly lower ve- 
locity dispersions of ~ 0.1 cm s -1 , satellite excitation of 
particle eccentricities in the last ~ 1 km of the ring is required. 
The presence of the shepherd satellites changes the precession 
condition for streamlines near the edges as well as maintaining 
the ring edge. Because shepherd satellites and ring edges are 
linked, the shepherding process has the potential to explain 
both the equilibrium ring masses and the ring surface density 
radial profiles. More work will be required before we under- 
stand the boundary conditions for a narrow, eccentric planetary 
ring. Whether such a ring can have a negative eccentricity gra- 
dient is presently unknown. We conclude that a modified CG 
model holds the best promise of finally making it possible to 
quantitatively reconcile observational mass estimates for nar- 
row, eccentric rings based on radio occultation constraints of 
ring particle sizes with the dynamical mass estimate based on 
the observed ring state. 
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